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Abstract 



We have calculated the 3-point correlation function for the elec- 
tromagnetic interaction of the pion in an ensembles of instantons and 
anti-instantons, modelling the QCD vacuum. The results are well de- 
scribed by a pion pole and the pion formfactor is extracted, nicely fol- 
lowing a standard monopole fit. The experimental data on the formfac- 
tor are well reproduced, provided an average instanton size is fixed at 
p = 0.35 ± 0.03 fm, the same as found for a variety of other correlation 
functions and was found on the lattice. 
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I. INTRODUCTION 



Significant progress has been made during the last few years in understanding of 
the QCD vacuum and hadronic structure in terms of instantons. In a series of papers 
it has been shown that a specific model, which assumes ensembles of instantons 
and anti-instantons in the QCD vacuum, can not only describe the gross features of 
the vacuum, such as spontaneous chiral symmetry breaking or topological susceptibil- 
ities)], but can also quantitatively describe a large number of mesonic and baryonic 
correlation functions, in agreement with phenomenlogy || and lattice calculations 0. 
Furthermore, the parameters of the ensemble have been confirmed by lattice studies 
an d also the dominant role of instantons was directly demonstrated by removal 
of all other kinds of the gauge fields from the configurations by the so called cooling 
procedure 0. 

In view of those developments, wider practical applications of the instanton-based 
models are justified. Since studies of 2-point correlators have produced hadronic masses 
and coupling constants in agreement with data, it is natural to perform more detailed 
studies and test whether those models can or cannot describe details of hadronic struc- 
ture as well. Three point correlation functions can describe hadronic couplings to ex- 
ternal fields, and for this aim they have been considered within the QCD sum rule 



approach (e.g. 0), the light-cone formalism (e.g. pjlJTTI ), the generalized impulse 



approximation within non-perturbative Dyson-Schwinger equations []T2 1 and of course 
on the lattice JT3HT5J. 

Only very recently the first (and the simplest) 3-point function related with the 
pion formfactor has been studied in the so called single-instanton approach by Forkel 
and Nielsen ||16|| . For Q 2 ~ 1 GeV 2 they have found good agreement with data. 



1 See e.g. 0] for a recent review of chiral symmetry breaking mechanism by instantons 
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Significant advantage of their approach over the previous QCD sum rule calculation 
is that the pseudoscalar (rather than the axial) currents are used in order to create 
pions. It makes a significant practical difference, because the pion is coupled to the 
former current with large constant and to the latter with the small coupling f n . 

Technically Forkel and Nielsen follow closely ref. [IT] in which the 2-point pion 
correlator was discussed in the single-instanton approximation: in it the propagators 
are calculated using the zero-mode solution of the Dirac operator in the presence of 
one close instanton (or anti-instanton), whereas the effect of the remote instantons is 
taken into account in some mean field approximation. This is a good approximation, 
provided all relevant distances (x, y) are short compared to inter-instanton separation 
R ~ 1 fm >> (x,y) ~ 0.2 fm. In practice, it is a very important limitation, since 
it does not allow to separate the pole term from the non-pole or "continuum" con- 
tributions unambigously. Therefore Forkel and Nielsen had to extract the formfactor 
by a complicated fitting procedure, similar to what is done in the context of QCD 
sum rules. Also, they have performed a traditional Borel transform of the correlator, 
while in fact the whole analysis can better be made directly in (Euclidean) space-time 
represent aion, in which all relevant formulae are much simpler. 

In the present paper we also study with better accuracy a variety of two point 
correlation functions, and check how the obtained masses of hadrons are correlated 
with hadronic formfactors. The electromagnetic form factor of the pion is of special 
interest, because the pion is the Goldstone boson of the spontaneously broken chiral 
symmetry and therefore the quark attraction in this channel is very large. Furthermore 
in the models used the pion is bound by instanton-induced forces, and in a way its 
formfactor is basically the instanton formfactor. 

We are using the methods developed in |] which allow to evaluate quark prop- 
agators in a multi-instanton background numerically, up to rather large distances. 
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Therefore, we are not restricted to small x, y << R and thus have no significant 
"background" from non-pion states. As a result, we can extract the pion formfactor 
with much better accuracy, and even study how it depends on various modifications of 
the underlying ensemble. Since the pion formfactor is also relatively accurately mea- 
sured, those studies proved to provide a very rigid constrained on such microscopic 
parameters such as the mean instanton size. 

We use two formulations of the statistical instanton-based models: the simplest 
Random Instanton Liquid Model (RILM) |T8|,p]fl] and the much more developed Inter- 
acting Instanton Liquid Model (IILM) (TJ5],|2(J. While the former one has simply the 
previously determined total instanton density and average size as input and assumes 
a random distribution in all collective coordinates, the latter one is a statistical model 
which predicts correlations between the instantons in the liquid. Some of those corre- 
lations are known to be qualitatively important: e.g. quark-induced interactions lead 
to a screening of the topological charge and correct the large distance behaviour of the 



correlators in channels, where the random model produces repulsion |21]| . It is also 
important that in the IILM the instanton parameters were not fixed but determined 
selfconsistently, and after the instanton interaction is fixed the only parameter left is 
Aqcd- If its value is fixed from an average density of n = 1 fm~ , an average instanton 
size which emerges is somewhat larger than the original p = 1/3 fm [113 or recent 
lattice results 

The organization of the paper is as follows. In section II. we summarize the phe- 
nomenlogical knowledge of the 3 point correlation function for the pion. In section III. 
we show how to derive this correlator within the Instanton Liquid model and show also 
the various short distance approximations to the full simulation. Sect. IV discusses 
the numerical procedure and results for the RILM and sect. V for the IILM. In sect. 
VI we make some concluding remarks. The App. A. gives some formulas for the direct 
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instanton approach and App. B gives formulas for the continuum contribution of the 
3 point function. 



II. ELECTROMAGNETIC PION FORM FACTOR - NOTATIONS AND 

PHENOMENOLOGY 

The electromagnetic structure of the pion can be measured in the time-like region, 
the p-meson pole containing part of the formfactor or structure function and in the 
space-like part. In both cases the formfactor depend only on the four- momentum 
transfer of the pion which we denote by Q 2 = —q 2 . (Below we always use space- 
like momenta, Q 2 > 0.) While the time-like region can be investigated in collider 
experiments such as e + e~ — > 7r + 7r~, the space-like region is obtained from scattering of 
pions off atomic electrons (for small q 2 ) or via the production of pions in e~p — > e~n + n 
(for larger q 2 ). 

From the theoretical point of view the interaction of the hadrons with external 
currents like the electromagetic one is most conveniently written in terms of three 
point correlation functions. The three points denote the creation and annihilation 
points of the hadrons as well as the location of their interaction with the external 
current j®(y)' 

n M (p, q)= J d'x J A e*" e^(0 | TjF{-x/2)tf(y)ft{x/2) | 0) (1) 

where p = (p + p')/2 is the average of the momenta p,p' connected to the pseudoscalar 
currents. The charged pions are represented by the quark bilinears 

/ \ 



d{x)i^u{x) 



u(x)i-y 5 d(x) j 



-^q{x)i^T ± q(x) := g(x)r>g(x) (2) 



and jfjfty) = q(x)Q-/ fl q(x) = Ei=u,d e i q i (y)') '^(y) is the electromagnetic current. 



The 3-point function in the r.h.s. of eq. (p]) is the quantity we will calculate below. 
Its transform, the l.h.s. of eq. flU) , can be related to physical observables by a double 
dispersion relation p2| p| 



over the spectral density p ll (s 1 , s 2 , Q 2 = -q 2 ) = pi(si, s 2 , Q 2 )p/i + P2(si, s 2 , Q 2 )q^,- As 
explained below, the e.m. formfactor of the pion can be extracted from the former 
structure containing p M and pi(si, s 2 , Q 2 ). 

A Borel transformation in p and p' was used in earlier papers, with the purpose 
to eliminate the polynomial subtraction terms. Our philosophy is however different: 
we make an inverse Fourier transfer to coordinate representation instead, in which the 
propagators in the instanton background are naturally obtained. As shown in details in 
H for the 2-point functions, the coordinate space is as good as Borel representation for 
all our purposes.^ The second reason for using the Borel transform in QCD sum rules, 
namely to suppress the larger mass resonances, is also not needed here, because in the 
coordinate-space correlator these states are automatically suppressed at distances we 
will work with. 

The spectral density p M (si, s 2 , q 2 ) is proportional to the imaginary part of the po- 
larization operator. We approximate it as the pion contribution given by 

2 The double integral is needed because one has to use a complete set of physical states 
both in the incoming and the outgoing channel. For a recent discussion on the dispersion 
relations cf. ^J24j|. Although the eq. (|3[) is not well defined due to the divergent behaviour 
of the dispersion integral for large momenta s,s', one can improve that by subtraction of 
polynomials functions oip 2 ,p' 2 ,Q 2 . 

3 As we show in App. B., the coordinate space representation of eq. (^) is not only finite 
but reproduces also the correct short distant behaviour. 




(3) 
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P P oie{si )S2l q 2 ) = (27T) 3 5(s 1 -ml)5(s 2 -ml)(0 | j 5 (0) [ 7r><vr | jfj | tt)(0 | j 5 (0) | tt) (4) 
Here the overlap matrix element 

(0 | h{x) | nip)) = K(27r)-V 2 e^ (5) 



is given in terms of A„-, which may be represented as |]T9 



with f n corresponding to the experimental f n = 93 MeV. The definition of the elec- 
tromagnetic form factor F n (q 2 ) eq. (f|) is 

(7r(p) I j M (0) | Tr(p')) = eMq 2 ^, (7) 

which is normalized by F n i0) = 1. Then the eq. (^) can be written as 

n "fo q) = hp - q /2f + ^mp + q /2y + ^ + n — fe^) ( g ) 

with free quark continuum H^contipi l) ( c f- App. B), starting at some threshold value 
The experimental measured form factor in the space like region has been shown to 



be well parametrized by a monopole form factor [25,26], resembling is some way the 



exchange of a p meson in vector dominance models. The parameter in this fit with 



my = 679 ± 19 MeV |25j is however below the p meson massg The measurements are 
performed over the large range R of space-like Q 2 = 0.18 — 9.77 GeV 2 . In coordinate 



4 There is also a precise measurement of the formfactor for smaller Q 2 < 0.30 GeV 2 |p6[ , 
which yields, also from a monopole fit, my = 736 ± 9 MeV and is still below the p mass. 

5 See also the CEBAF proposal [^7]] for a future high statistics measurement of the formfactor 
via electroproduction at Q 2 ~ 0.5 — 5 GeV 2 . 
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space this measures the region around 0.05 — 0.5 fm. The upper limit of which is clearly 
dominated by non-perturbative effects in the QCD vacuum, as will be shown, and in 
general cannot be described by the operator product expansion. 

The fourier transform of eq. (H) is then conveniently calculated from 



II , , / A I 1 A ? I V - X ( - 5 - a ) I) / 2 , 2 ft \ 

ii UiVO i e {x,y) = —i I dq I da : ; — \ mi + q z a(l — a) 

M,p v ,yj Jo H Jo 32tt 4 I y - x(0.5 - a) \ V * H V ; 

2A^ e7r x. 



\ v / 1 + q 2 /rriy \ x \ 

where J\,K\ are the normal and modified bessel functions. The rather cumbersome 
formula for the continuum is given in App. B. In eq. @j we used already the monopole 
parametrization of the form factor [p5| |, which was shown to be a good description in 



the space like region of the form factor. 

III. THE INSTANTON LIQUID 

The 3-point function in question can be evaluated using the same quark propagators 
in the multi-instanton background fields as used for the 2-point functions in [ p!9| , p0[ 



and we will not repeat any details here. The generating functional of the IILM can be 
written as 

1 f N+ ^- -/(frf)/ f)(i)fp + m f -sT)(i )2 /)l"\w 



1 



f 

We will only shortly explain the meaning of the various expressions. The D is the 
fermionic Dirac operator in the presence of instantons and anti-instantons. These in- 
stanton solutions are described by a set of collective coordinates which are the color 
orientation Ui, the position Z{ and the size pi of the instanton. Therefore the original 
gauge measure T>A^ becomes essentially an integral over the collective coordinates dVti. 
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The d(pi) in eq. (|T0|) is the semiclassical instanton amplitude, which was originally 
calculated by 't Hooft: its two-loop version is used in [|i9| . 

The rrif are the explicit (chiral and vector) symmetry breaking current masses 



and Si n t denotes the classical interaction of the instantons, which is as usual |19 
approximated by a two-body interaction. 

In the case of the electromagnetic formfactor one has = / y fl Q a T a , where Q is 
the electromagnetic charge matrix. The spin-isospin matrices T^, (0 = ir,8), in eq. (0) 
correspond to rysT 1 * 1 / v^2 for the charged pions 7r ± and to r ± / y/2 for the isovector scalar 
5-meson. The reason to discuss also the 5 shortly will become clear in the following. 
From the generating functional eq. (|H]) a three point correlator for a mesonic state is 
obtained from 

TT ( \ — (V \ af5 (V V s — ^. £ 4 Z[t7,77, s°] 

" {X,V) " 1 * )ab{ ' )cd 5s^y)S m (-x/2)S^(-x/2) SffS{x/2)S4(x/2) ^ 

where latin indices correspond to spin indices and greek symbols to isospin indices. 

Since all currents of interests are color singlets, color indices are implicitly understood. 

In terms of the quark propagators the general form for the 3-point correlator, using 

isospin symmetry for the up and down quark propagator S u = Sd, follows as 

n M (x, y) = (e u - e d ) «trS(-x/2, y)^S(y, x/2) l5 S(x/2, -x/2) l5 )- 
(txS(y, y) ltl tvS(-x/2, x/2) l5 S(x/2, -x/2) l5 ) - 
{trS(y, y)7^r5(-x/2, -x/2) lh tr5(x/2, x/2) 75 » (12) 

Since the quark propagators in eq. (|1^) are evaluated in a given instanton anti- 
instanton background, averaging over different instanton ensembles is implicitly under- 
stood by the brackets in eq. (|T2 ) . Here the first term on the RHS of eq. ( |T2"D is the so 



called connected contribution and the second and third are the disconnected diagrams. 
One should note of course that the quark propagators are the full propagators in the 
background of the multi instanton/anti- instanton configurations so that the expression 



'disconnected' refers only to the 'quark lines'. One may note that the disconnected 
terms are in general obtained from the fermion determinant of the generating func- 
tional eq. fllOD Q . 



In the present approach there are two things to mention. In the coordinate repre- 
sentation and using the free, massless quark correlator]] 



/■yi ry ry ry rt 

^-U) = ^\ru(--)u ( -)\ ) = --± t (is) 



this second and third terms in eq. ( O ) would have a short- distance singularity. In the 
case of the vector current, one has to subtract the free propagator The discon- 

nected parts can be shown to vanish (which follows already from current conservation 



[29 1) and only the triangular diagram of eq. (|T|) survives. 

For very small Euclidean space-time separations | x |, | y |<C 0.2 fm the 3-point 
function H^x, y) is governed by the free Dirac propagator eq. ([13]) . For x ■ y — 0, 
which we use, it is especially simple 

lim U u (x,y) = —?r—iri—^ n . .„ (14) 

x, y ^o M ,yj 2tt 6 a; 4 (y 2 + x 2 /A) 3 v ; 



6 In the case of effective quark or semi-bosonized quark theories like the NJL model [28 
these expressions are known as the UV divergent one loop contributions, representing the 
polarization of the vacuum. 



7 Note that we are using anti-hermitian 7-matrices, so that eq. (13) is indeed real in Eu- 
clidean space. 

8 In the presence of the instanton also the non-zero modes get distorted [3C,|] and some 
disconnected contributions still survive and are even divergent Q. They however disapper 
after insertion of a path-ordered exponential q(y)Qj ll P exp (— i J^ +e A fJ/ (z)dz fJ/ )q(y + e) which 
makes the electromagnetic current gauge invariant and conserved. 
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For larger distances the non-perturbative effects come in. In the leading order for small 
distances the vacuum dominance approximation would suggest the following correction 
to eq. flU) 

5 (-f.f) = -lT (15) 

It is clear that the term of the order O (| uu |) vanishes due to spin traces, while the 
0(\ uu | 2 ) contribution in this approach reads^] 

y) = -^ff l (««> l 2 L + I w + ^) ( 16 ) 



Compared to eq. (fbl| ) , this next to leading order term has the same sign and indicates 
therefore an enhanced signal in this channel. This is similiar to the pion two-point 
correlator itself || and can be traced back to the attractive forces inside the pion. 

However such "vacuum dominance" approximation works only qualitatively and for 
small distances only. To go beyond those one needs a model for the non-perturbative 
effects. The basic feature of the instanton-based models is related to the zero mode 
solutions of the Euclidean Dirac operator D, which are obtained from 

D4>x = A^a, (17) 

for A = and exist, if the topology of the gauge field configurations is non-trivial. 
This is the case for the instanton solutions. From these solutions we approximate the 
full Dirac propagator S(—x/2, x/2) in the instanton liquid model by a sum of the zero 
mode contributions |],||] Szm(—x/2,x/2): 

W-f , x - ) = E<M-f )(/ 1 I J>ri(f ) (is) 

* I j j I D + m * 



and the non-zero mode term Snzm(—x/2,x/2) [ j5D| , pT 



3 Note that it is different from the 0(a s (^>^/) ) contribution in standard OPE analysis ||. 
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s 



NZM\ 



X X, 

2' 2' 



'2' 2' 



E 



SA- 



2' 2' 



5nf- 



X X, 

2' 2' 



(19) 



where Sj(—x/2,x/2) is the non-zero mode contribution in the presence of a single 
instanton^ 

Before discussing the result of the full calculation with eq. (|18| ) and eq. ( |19|) in 
eq. ( |12"D in an ensemble of instantons, it is instructive to investigate the effect of a 
single instanton, as done by Forkel and Nielsen. Then the general form of the zero 
mode propagator is given by 



S 



x x * A (_D*t(f; 

A 



(20) 



A + im 
where * A (a?) = E/LiC/*o 

(re — zi) can be expanded in terms of the zero modes ^o( x ) 

of individual instantons. Now one can argue |32| that for small distances the instanton 

i* closest to I and — | dominates so that the propagator can be expressed in terms of 

this single zero mode for an instanton/anti-instanton as 

c x ,2 



2' 2' 



E 



A + im 



y/ 



8m 



\ 



r + r~P R 



:-f-*>#(f-*> 



(21) 



where a dynamically generated effective mass m can be defined from eq. (pop by 

C/ I 2 



1 

m 



E 



A + 



(22) 



and 



P 



71 \ x \ (x 2 + p 2 ) 



,2\3/2- 



(23) 



10 This expression eq. (|TJ|) should be understood as the leading term in a multiple scattering 
expansion, cf. § for details. 
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In this approach the effect of distant instantons and chiral symmetry breaking is there- 
fore effectively included by replacement of the current quark mass m by fh, the dynam- 
ically generated effective mass |l7j], representing chiral symmetry breaking by other 



instantons. This was done by Forkel and Nielsen |L6| who has used the correlator 
y) with eq. fl2~l~|) . In the short distance limitQ this reduces to 

lim UJx,v) = — -^7 — — — — %>1 n . NO (24) 
*,i/-o M ,yj 5vr 4 fh 2 p 4 {y 2 + x 2 /A) 2 K ' 

Note that using the mean field results [0,0 for the constituent quark mass R and the 
condensate, namely (uu) ~ y/E/ p and fh ~ V^p, the eq. (E4p gives a correlator which is 
parametrically larger than eq. ( |i~6|) by the inverse packing fraction 1// ~ l/(np 4 ) 3> 1. 

However for the full Instanton Liquid simulation the formulas of Forkel and Nielsen 
are consistent to our approach for short distances. But neither the vacuum dominance 
nor the single-instanton approach |RJ can describe the correlator for larger distances 
x,y ^> 0.2 fm necessary for a clear separation of the pion contribution. This goal will 
be reached in the next section, by using the numerically calculated correlators. 



IV. CORRELATORS IN RANDOM INSTANTON VACUUM 

Numerically the propagators are calculated as solutions to eqs. (|T8|), (|T9"D. Further- 
more for the RILM we took 256 randomly placed and oriented instantons (half of them 
anti-instantons) into a periodic box of (5.67) 2 x (2.82) 2 fm 4 , where the long box lengths 
correspond to the x, y directions, in which the correlator is actually measured. The 
numbers are chosen to reproduce an instanton density of h = 1 fm -4 . Furthermore, 



See our discussion of the full expression in the next section and App.A. 
12 Note that it is not equal to the so called constituent quark mass (and it is in fact about 
twice smaller) because the effective mass in question is not evaluated at zero momentum. 
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we take a variable average instanton size ( p = 0.28 fm, 0.35 fm, 0.42 fm) in order to 
show the sensitivity of the form factor to the instanton size. We have evaluated 2- and 
3-point functions and studied them in detail. 

In Fig.l we show the two point correlation function. According to our pole plus 
resonance Ansatz (eqs. |A5| , |A6| ) for the correlator we give in the figure the best fit 
for the pion mass m n and the coupling constant X n . Because of our finite box we 
used a current quark mass of m u = = 20 MeV which is somewhat larger than the 
physical quark masses. Using the Gell-Mann Oakes Renner relation, the measured pion 
mass values in Fig. 1 are scaled down to the physical pion mass, see Tab. (|) . The 
extrapolated value is close to the experimental value. For the same configurations we 
determined then these parameters by fitting the three-point correlators shown in Fig. 
2. In principle the correlator is a function of two points, x and y, which we have chosen 
to be orthogonal. So one can separately check that the form factor has no or neglible 
x dependence and then look for the y dependence. The x value has to be chosen large 
enough, so that the correlator is clearly dominated by the pion. However is x too 
large the statistical errrors become too large and no extraction of the parameters of 
the form factor is possible. Therefore in practice we followed an alternative scheme, 
which covers both situations simultaneously and furthermore saves some computing 
time: we fixed the ratio of x and y to x = 2y. In app. B we also give formulas for the 
continuum contribution, though we want to stress that our parameter fixing is totally 
independent of this parametrization since we are considering distances x ^> 0.5 fm, 
where these contributions are highly suppressed compared to the pole contribution. 

The pion parameters, m n and X n , found from 2- and 3-point correlators agree within 
the error bars. Established consistency, we have further determined the mass in the 
monopole formfactor, which actually can be done rather accurately. We have found 
that this mass (or the pion size) turns out to be sensitive to the instanton size, see 
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Fig.3. Remarkably enough, the experimental corridor (two horizontally dashed lines) 
indicate the preferred instanton size to be p = 0.35 fm, in approximate agreement with 
1/3 fm |L8| and in very good agreement with lattice results [0 pointing out p = 0.35 fm 



as well. 

As can be seen in the figure a further improvement of the experimental measured 
formfactor, as it is planned [P^j , and increasing the accuracy of the instanton cal- 
culation could provide a very powerful tool actually to measure the instanton size. 
This is a remarkable statement since it directly connects a physical observable like the 
formfactor to an intrinsic property of the QCD vacuum, the size of the instantons. 

In addition in Fig. 7 we compared our full RILM result to the one-instanton ap- 



proach of the pion formfactor of Forkel and Nielsen |L6| , but now plotted in coordinate 
space and without Borel transform. As can be seen, this approach is reliable for 
distances not larger than 0.4 — 0.5 fm. The same qualitative behaviour of the two 
approaches was found for the pion correlator itself ||34j| . 

V. CORRELATORS IN INTERACTING INSTANTON VACUUM 

On general grounds, the pion channel is a strongly attractive one, and it was 
expected that pion properties (including form factors) do not depend on the details of 
the instanton ensemble, such as correlations etc. 

In order to check these expectations we repeated the calculations for the IILM. Our 
results are shown in Fig. 5 and 6 for 2- and 3-point correlators (analogous to Fig. 1,2). 
Because of the particular interaction assumed, the ensemble has an average instanton 
size of p = 0.42 fm |y|, so we are not able to compare different sizes for this model. 



However, by comparing the IILM with the RILM for this common size of p = 0.42 fm, 
we can find out the role of the correlations. As can be seen from Tab. (||) , we have in 
fact found that the pion coupling constant depends rather strongly on the model (to 
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a degree, this may be traced back simply to different values of the quark condensate), 
while the pion and the monopole mass (of the form factor) are less dependent. Within 
the error bars of our fits, the m* and my of Tab. ([H]) seem to agree with the values 
of Tab. (|) for p = 0.42 fm. 

The same can be done for the isovector scalar meson, the S, as shown in Fig. 
4. However, as already mentioned, in this channel the repulsion is very large and the 
simple random instanton liquid does no longer give a good description of the correlator. 
On very general grounds, the correlator should be positive for larger distances, whereas 
the RILM correlator changes sign for x = 2y ~ 0.6 fm and crosses the axis again 
at x = 2y ~ 1.4 fm. The situation changes dramatically if we consider the IILM 
correlator. This one, as can be seen also in Fig. 4, stays positive over the whole range 
of available distances. However for x = 2y ~ 1.8 fm the error bars become larger 
and the correlator increases again, which is an unphysical effect. Since the mass of 
the 5 is too large to be measured precise enough in our model, we do not attempt to 
determine the mass of the possible form factor. The figure should serve only as a further 
qualitative justification of our IILM model, in which the correlation and interaction of 
the instantons provide the right behaviour for the strong repulsive channels. 

The last issue we address in this work deal with the old question of vector domi- 
nance. We remind the reader that it suggests that a complete pion formfactor is given 
by a rho-meson pole alone, and thus the fitted mass in the formfactor my is nothing 
else but m p . It is certainly approximately true in nature: and therefore we have tried 
to check whether indeed there exist a strong correlation between them (in the models 
used). Therefore we have determined also the parameters of the p meson correlator 
(shown in Fig. 8). However we found that: (i) m p is only very weakly dependent on 
the instanton size, from 875 to 915 MeV in our given interval for the instanton sizes; 
(ii) m p is more sensitive to correlations than my; and (iii) my is also significantly 
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(more than 10%) smaller than the experimental m p We therefore conclude, that 
in the instanton-based models the mass in the pion form factor my has nothing to do 
with the p-meson mass, and the fact that they are close numerically is probably just 
a coincidence. 



VI. SUMMARY AND DISCUSSION 

We calculated the pion electromagnetic formfactor for space like momentum trans- 
fer, based on a three-point Euclidean correlation function of two pseudoscalar, isovector 
currents and an external electromagnetic current. Two variants of the instanton-based 
QCD vacuum models were used, a random one, RILM, and the interacting one, IILM. 
Our main result is the existence of the direct connection between the size of the in- 
stanton and the size of the pion. 

In contrast to earlier work ||16|| , we have calculated the correlation function in 



coordinate space at sufficiently large distances, clearly separating the pion pole con- 
tribution from those of non-resonance states. We have accurately determined the pion 
coupling constant A w , the pion mass m w and the parameter my of the monopole fit to 
the electromagnetic formfactor. The first two parameters are fitted both from 2-point 
and 3-point correlators, and the results are consistent within errors. 

We have calculated the monopole mass of the formfactor my for three different sets 
of average instanton sizes at fixed instanton density n — 1 fm~ 4 and find that it is in- 
deed sensitive. Thus the pion size is directly related to the instanton size. Furthermore, 
the experimental monopole mass my ~ 679 ± 18 MeV is only reproduced for a mean 
instanton size around 0.35 fm, a value consistent with direct lattice measurements |7]]. 



13 The latter one is in the instanton model unfortunately too large, probably due to the 
crudeness of our phenomenlogical Ansatz. 
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As a further check we used the streamline version |2l]{ni of the interacting instanton 



model (IILM). In this case the inverse monopole mass is small 0.16 fm, or my ~ 
1200 MeV. It does not agree with data but agrees well with the RILM result for the 
same mean instanton size of p = 0.42 fm. These results indicate that the formfactor 
does not depend very much on the correlations between instantons. Other quantities 
like the condensates, e.g., do however depend rather strongly on which model is used 
and therefore on such correlations. So, the IILM can and should be improved. 

We have considered the issue of vector dominance, suggesting that the formfactor 
mass my is nothing but m p . However, in our model used, we have found that both 
quantities show different dependence on the instanton size and correlations in the 
ensemble. We therefore suggest that vector dominance should not hold in general and 
seems to be more a kind of a coincidence. This will be checked by considering other 
3-point functions in forthcoming publications. 
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TABLES 





2-pc 
\l /2 [ MeV] 


int 

m* [ MeV] 


3-po 
\l /2 [ MeV] 


int 

m*[ MeV] 


m v [ MeV] 


(uu^^l MeV] 


[ MeV] 


p = 0.28 
p = 0.35 
p = 0.42 


455±5 
490±5 
520±5 


265±5 
285±5 
255±5 


455± 5 
480± 5 
510± 5 


265±5 
290±5 
265±5 


1250±50 

680±20 

550±20 


280.3± 0.46 
259.1± 0.25 
245.2± 0.39 


139 
150 
136 



U[ MeV] 



149.9±5 



106.6±5 



79.9±5 



TABLE I. The pion properties in the RILM for different values of the instanton size 
p. The pion coupling constant and pion mass m* for m u = 0.1 A ~ 20 MeV obtained 
from the 2-point and 3-point function are determined independently and coincide within the 

1/2 

uncertainty of the fit. For comparison, Forkel and Nielsson used \^ = 363 MeV, which 
either corresponds to a rather low value for the condensate (nn) 1 / 3 ~ 200 MeV or a large 
f n ~ 180 MeV. The physical pion mass m n in the table follows from an assumed Gell-Mann 
Oakes Renner scaling relation with m u + = 11 MeV. 





2-pc 
\l /2 [ MeV] 


int 

m*[ MeV] 


3-po 
\l /2 [ MeV] 


int 

m*[ MeV] 


m v [ MeV] 


(un) 1 / 3 ! MeV] 


ma [ MeV] 


p = 0.42 


395±5 


260±5 


370± 5 


295±5 


1200±50 


217.3 ± 0.55 


113 



U[ MeV] 



93 ±5 



TABLE II. The pion properties in the IILM for A = 306 MeV. The pion coupling 
constant A,r and pion mass m* for m u = 0.1A obtained from the 2-point and 3-point function 
are determined independently. The physical pion mass m w in the table follows from an 
assumed Gell-Mann Oakes Renner scaling relation with m u + = 11 MeV. 
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APPENDIX A: ONE INSTANTON FORMULAS 



1. The pion correlator 

The pion correlator for pseudoscalar currents j${x) for short distances in the asymp- 
totic regime is simply given by 

< | Tj 5 (x)j 5 (0) | >= ^-i := n(x) /ree (Al) 

Corrections to this correlator at larger distances are given by the direct instanton 
contribution: 

< | Tj 5 (x)j 5 (0) | > = -,n—n A 



7T 2 m 2 ^/Ax 2 p 2 + x 4? 
arctanh^= X 2 4 8 (x 4 + 2p 2 + 2p 2 x 2 ) + ^Ax 2 p 2 + x 4 -^ (2x 4 - 10x 2 p 2 - Yip 

■= U(x)inst (A2) 



where m (cf. eq. |22|) is the effective constituent quark mass |35| , |36| The asymptotics 



of the one-instanton formula eq. ( |A2|) for short distances is given by 

In 2 

limll(x)i ns t = -^—^p 7T ( A3 ) 
x^o ix 2 m 2 5p 8 

which is suppressed by the free contribution eq. ( [Al] ) with the power x 6 , so that the 
sum of Yi{x)inst and U(x) f ree , normalized with H(x) f ree , goes to unity. This is also the 
case for our instanton liquid ensemble of Fig.l and 5. For large distances the limit of 
eq. ( |A2| ) becomes 

12 1 2 

lim Yi(x) inst = —n—p 4 —— (A4) 
x ^oo tt 2 m 2 3p 2 x b 

Here the zero mode contribution is proportional to the free part eq. (|A1|) and actually 
dominates the sum of both expression. For completeness the pole contribution of the 
2-point function in a simple pole plus continuum Ansatz is || given by 
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A 2 1 

n(x) poie = —^—m^ | x | #1(771* I x |) (A5) 
Air z x A 

whereas the continuum results || reads 

n(x) con4 = (^(^(lGEx + A(Exf) + # (£x)(8(£x) 2 + (Ex) 4 )) (A6) 
Indeed the small size limit of eq. (|A6| ) for a fixed threshold E is given by 



Km Il(x) cont = Il(x) f ree (A7) 



whereas the normalized pole contribution, i.e. linx^o n(x) po i e /U(x) f ree , vanishes. 



2. The pion EM-formfactor 



From our general expression eq. ( |12"D and the one-instanton Ansatz eq. (^T|) the 
correlator is given by 

rr I \- ~ ie * mp4 1 ( > 

ii^x, V )- ^ _^ (2/2 + x2/4)2 y az 

y^x 2 /4: + 2yz - z 2 ) + x M /2(-y 2 + 2yz - z 2 ) + ^(-y 2 - x 2 /4) 
\y-z\((y-z) 2 + p 2 ) l - 5 ((x/2 - zf + p 2 ) 3 | z/2 + z | ((x/2 + zf + p 2 ) 1 - 5 1 j 

where the trivial averaging over color orientations is already done. The averaging over 
the instanton position in eq. (|A8|) can be performed after introducing 5 Feynman 
parameter integrals for the 5 denominators in eq. ( |A~8| ) . After this, it is even possible 
to perform 2 of the Feynman integrals analytically, which results in the rather awkward 
expression 

u-p{x,y) = ——5760-^- — fj— d a x \ d a 2 aa 3 a u / — 

7r 8 m z (y z + x 1 /Af Jo Jo Jo y 03 

tt 3 (1 - ai - a 2 ) 356 3 + 120a6 2 + 1446a 2 + 64a 3 
3840 a 4 - 5 (a + 6) 3 - 5 

(x 2 /4((ai - a 2 - a 3 ) 2 + (ai - a 2 - a 3 )) + y 2 ((ai + a 2 + a 3 ) 2 + (al — a 2 — a 3 ))) 

7r 3 (l - ai - a 2 ) 56 2 + 126a + 8a 2 " 



960 a 3 - 5 (a + 6) 2 " 



(A9) 
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where 

a = p 2 (l - a 3 ) + x 2 /4((ai + a 2 + a 3 ) - (ai - a 2 - a 3 ) 2 ) + y 2 ((ai + a 2 + a 3 ) - (ai + a 2 + CJ3) 2 ) 

(A10) 

and b = — p 2 (l ). Similiar to the 2-point function one can consider the 

short distance behaviour of this function and finds 

lim T (x, y, 0) = - ^ n (All) 

Because our triangular diagram for the EM form factor in this approximation actually 
consists of 2 zero mode propagators and one free propagator, the zero mode part 
reduces to a constant, similiar to eq. (|A3|) , and a single quark propagator in eq. ( [All ) 
remains. The full behaviour of eq. ( |A"9| ) can be seen in Fig. 7 compared to the full 
(RILM) result. Obviously the one instanton formula deviates from the full result for 
distances larger than 0.4 fm. 



APPENDIX B: CONTINUUM CONTRIBUTION 

For the continuum contribution one has to calculate the double discontinuity of 
the spectral density from the expression for the triangular diagram according to the 
Cutkovsky rule as 

N 1 r rl A k 

P»(sn s 2 , (?) = Ll—J —6(k )6(k 2 )6(k - Po )6((k - p) 2 )6(k - p' )8((k - p'f) 



(Bl) 



with the result 



PnK.su s 2, -q 



sis 2 N c 



P»q 2 -q»{s 2 - Sl )/A (B2) 



7T 2 A 1 - 5 

where A = (si + s 2 + g 2 ) 2 — As\s 2 . Inserting this into the dispersion relation and 
transforming the resulting expression into Euclidean space, one can use the Feynman 
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parameters as before for the pole contribution to perform the p-integration analytically. 
In order to model the continuum one usually introduces a threshold parameter sq into 
dispersion relations. If one does it, the continuum contribution reduces to a four 
dimensional integral, which can be easily done numerically. It reads 

N c Sl s 2 q 2 



rOC rl rOO 

n^, co „i(x,y) = -iN c / dsids 2 da dq (1 - 9(s 1 - s )6(s 2 - s )) 
Jo Jo Jo 

fM^-K^ x | M) 



64tt 6 A l 



, x 
s' — s 



Ji(q 


y + x(0.5 — a 


)l) 




y + x(0.5 — a) 





-2Ji(g j/ + x 0.5 - a) , , r, 

| y + x(U.o — a) | 



(B3) 



where 

M 2 = g 2 a(l - a) + sa + (1 - a)s > (B4) 



and A = (si + s 2 + q 2 ) 2 — 4siS 2 . In eq. ( |B3|) we introduced a simple model for the 
continuum by introducing the ^-functions [[24| . Below we argue why their detailed form 



is not important here. 

Now considering our preferred geometrical arrangement x y = it is clear from eq. 
( |B3| ) that the x^ part in the last two lines vanishes due to the a-integration. Finally 
the y^ terms vanishes because these are antisymmetric in s <-> s' 14 . Therefore only 
the first term in eq. flB"5|) finally contributes, so that 



n AtjPO i e (x, y) + Tl^ cont (x, y) = 

Ji(q | y - x(0.5 - a) |) 
| y — x(0.5 — a) \ 
2A 2 r e 7r x 



roc rl 

—if dq I da 
Jo Jo 



, q 2 

32tt 4 



m 2 + q 2 a(l — a) 



K\ ( | x | Jml + q 2 a(l — a) 



1 + q 2 /rriy \ x 



14 Using that M 2 is invariant under simultaneous change of s «-> s' and a <-» (1 
as the remainder of the formula eq. (|B3| ) . 
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— a) as well 



/"CO rOO 

+ / ds! / c/s 2 (1 - 6»(si - S )6 I (S2 - s )) 



7V C sis 2 ? 



(I x I M) 



2vr 2 A 1 - 5 

Ji (Q I y + #(0.5 — a) |) 



(B5) 



| y — x(0.5 — a) \ 

Similiar to the 2-point function eq. ( |A6| ) the small size limit of eq. ( [B"5| ) , which is 
actually the small size limit of n jUjCont (x, y), is given by H^x, y) of eq. (p^) . First one 
should note that eq. (|B5| ) is finite with neither referring to subtraction terms, which 
would be necessary in momentum space due to the insufficient fast decreasing spectral 
density, nor to Borel transformations. The latter method eliminates the subtraction 
terms, because these are known to be polynomials in either p or p' and therefore vanish 



after taking derivatives finite times |24| . In the present approach, the fourier transform 
of the dispersion relation into coordinate space, which is finite a priori, circumvents the 
introduction of any subtractions. However one should also stress, that the continuum 
contribution for a threshold of So > 1 GeV is highly suppressed for distances above 
1 fm. We have given those expressions for completeness and their inclusion does not 
affect our parameter fixing at all. 
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75.0 



pion 2-point function 

different instanton sizes (lambda.mpi) 



/ 




0.0 1.0 

x[fm] 

Fig. 1: The two-point correlation function of the pion, normalized to the free correlator 
and for three different instanton sizes in the random configuration (RILM). The instanton 
density is kept constant n = 1 fm -4 . The parameters of the pole contribution, and m^, 
are given in the legend. 



28 



pion 3-point function 

electromagnetic formfactor (MV,lambda,mpi) 
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Fig. 2: The three-point correlation function of the pion and an external electromagnetic 
current, normalized to the free correlator and for three different instanton sizes in the 
random configuration (RILM). The instanton density is kept constant n = 1 fm -4 . The 
parameters of the pole contribution, my, \ n and in this order, are given in the legend. 
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Fig. 3: The squared electromagnetic radius of the pion as a function of instanton size for 
the random ensemble (RILM). The instanton density is kept constant n = 1 fm -4 . 
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delta 3-point function 
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Fig. 4: The three-point correlation function of the delta (scalar, isovector) and an ex- 
ternal electromagnetic current, normalized to the free correlator in the random configura- 
tion (RILM) and the interacting ensemble (IILM). The instanton density is kept constant 
n = 1 fm~ 4 . 
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pion 2-point function 

(lambda.mpi) 

30.0 | . 1 . , 




x[fm] 

Fig. 5: The two-point correlation function of the pion, normalized to the free correlator 
and in the interacting configuration (IILM). The average instanton size is p = 0.42 fm and 
instanton density is p = 1 fm - . 
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Fig. 6: The three-point correlation function of the pion and an external electromagnetic 
current, normalized to the free correlator and for the interacting configuration (IILM). The 
averaged instanton size is p = 0.42 fm and the instanton density is p = 1 fm - . 
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Fig. 7: The three-point correlation function of the pion and an external electromagnetic 
current for the RILM and the one- instanton (linst) formula for three different instanton 
sizes. 
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rho 2-point function 

different instanton sizes (mrho.EO, lambda) 
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Fig. 8: The two-point correlation function of the p-meson for the RILM and three different 
instanton sizes, normalized to the two-point correlation function. In the legend, the p meson 
mass, the threshold parameter Eq and the p coupling constant are given. 
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